Competing Ferri- and Antiferromagnetic Phases in Geometrically Frustrated LuFe 2 04 
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We present a detailed study of magnetism in LuFe204, combining magnetization measurements 
with neutron and soft x-ray diffraction. The magnetic phase diagram in the vicinity of Tjv involves 
a metamagnetic transition separating an antiferro- and a ferrimagnetic phase. For both phases the 
spin structure is refined by neutron diffraction. Observed diffuse magnetic scattering far above Tjv 
is explained in terms of near degeneracy of the magnetic phases. 
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PACS numbers: 75.85.+t, 75.60.Ej, 75.25.-j, 75.30.Kz 

Magnetoelectric multiferroics are of interest for novel 
storage devices [l|, Q. LuFe204 was proposed to be a 
multiferroic with a novel mechanism for ferroelectric- 
ity, based on Fe 2+ /Fe 3+ charge order (CO) [3J below 
Tco'~ 320 K. Mainly for this reason, but also due to un- 
related effects such as giant coercivity, it is currently at- 
tracting a lot of attention Both charge and spin 
degrees of freedom are localized at the Fe sites, which 
are contained in triangular Fe-0 bilayers, a highly frus- 
trated arrangement. For the CO, competin g in stabilities 
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were in- 



suggested by diffuse scattering above T ( 
deed linked to geometrical frustration [6|. Similar geo- 
metrical frustration effects can also be expected for the 
magnetism, the elucidation of which is important for un- 
derstanding the magnetoelectric coupling and other in- 
triguing effects such as giant coercivity |la |. 

LuFe204 typically exhibits magnetic order or freezing 
below about 220 - 240 K. There is consensus that the 
Fe spins have a strong preference to be aligned || Chcx, 
perpendicular to the layers 1^-21|. The magnetic be- 
havior thus arises from Ising-spins on triangular lattices. 
Consistent with the highly frustrated arrangement many 
unusual effects have been observed in different samples, 
including various cluster or spin glass states 
magneto-structural transition at T^r~170K [R 
an anomalous "field- heating-effect" [20|. Strong sample- 
to-sample variations in magnetic behavior are found, at- 
tributed to tiny variations in oxygen stoichiometry. De- 
spite the high current interest, the details of the magnetic 
field (H) - temperature (T) phase diagram underlying 
these unusual behaviors have not yet been established. 

In this work, we present a detailed study of the H — T 
phase diagram above Tlt of LuFe2 04 [22J by magneti- 
zation and neutron and soft x-ray diffraction, revealing 



competing antiferromagnetic (AFM) and ferrimagnetic 
(fM) spin structures. The main focus is on samples with 
sharp magnetic transitions at Tjv ~240 K to long-range 
spin order [l5| , which we propose to best approximate the 
intrinsic defect-free magnetic behavior of LuFe204. We 
demonstrate that at T/v and H — fM and AFM instabil- 
ities, which correspond respectively to ferro and antiferro 
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FIG. 1. (Color online) (a) Magnetic field H - temperature T 
phase diagram, which exhibits a paramagnetic (PM) an an- 
tiferromagnetic (AFM) and a ferrimagnetic (fM) phase, ex- 
tracted from various M(H) and M(T) curves. The hysteretic 
region where either fM or AFM can be stabilized is hatched. 
Arrows across phase lines indicate for which measurement di- 
rection it is observed given the hysteresis. Spin structure in 
C2/m cell [H of the AFM (b) and fM phase (c). Grey arrows 
indicate bilayer net magnetization. 
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FIG. 2. (Color online) T-dependence of various properties, 
all measured on cooling, (a) Magnetization M, compared 
with another sample. Inset: Inverse susceptibility H / M with 
Curie- Weiss fit (dashed blue line) from 450 K to 750 K. (b) Ac- 
susceptibility measurement with different driving frequencies; 
the inset shows an enlarged area at Tjv. (c) Integrated in- 
tensity of (00|) x-ray reflection at 706.4 eV (Fe-I/3) and the 
(||0) neutron reflection, both in H = 0. Inset: Energy scans 
across the Fe-L edges at (00 5) with different incoming and 
outgoing polarization directions, (d) Integrated intensity for 
the (|||) neutron reflection and M(T), both in 2.5T. 



stacking of equivalently ordered bilayers, are nearly de- 
generate. These bear a striking resemblance with the 
two nearly degenerate CO instabilities [f| [l2[ at Tco, 
which we attribute to the similarity of binary (Ising spins 
or valence states) order emerging from competing inter- 
actions on the same strongly frustrated lattice. Diffuse 
magnetic scattering above Tjv indicates a random stack- 
ing of still individually ferrimagnetically ordered bilayers. 
We emphasize that although AFM-fM meta-magnetism 
has not been reported previously and may not be resolv- 
able in the majority of LuF^CU samples, our results have 
strong implications for the general nature of magnetism 
in this material. In particular, our results underline the 
importance of geometrical frustration in LUFC2O4, both 
for charge and spin order. 

We studied various LuFe2 04 single crystals from the 
same batch as in 



12Hl5|. Dc magnetization M and ac 



susceptibility \' measurements in if ||ch ox were performed 
with commercial (Quantum Design) equipment. Polar- 
ized neutron diffraction in if ~ (except a small guide 
field less than 10 Oe) was performed on DNS at FRM-II 
and non-polarized neutron diffraction with if ||ch cx up to 
2.5 T on D23 at ILL and C5 at Chalk River Laboratories, 
all using the crystal labeled S2 in 15j. Resonant x-ray 
diffraction at the Fe L3 edge was performed on the SIM 



beamline (RESOXS endstation) at the SLS. For compar- 
ison with previous work all reflections have been indexed 
in hexagonal notation. 

M(T) measurements show variations of magnetic prop- 
erties even among samples from one batch. One ex- 
treme exhibits characteristics matching those of [23j |. 
with a strongly frequency-dependent peak in \' around 
225 K indicating a transition into a glassy state. The 
other extreme exhibits at T/v~240K a sharp peak in M 
and x' without frequency-splitting (Fig(2^,/b) and shows 
(FigUa) sharp magnetic reflections in neutron diffrac- 
tion, indicating 3D long-range spin order rather than a 
glassy state. We also characterized samples with interme- 
diate properties: the peak signifying magnetic ordering 
is shifted to lower T (Fig. [2^, dotted line) and becomes 
weakly frequency-dependent. This indicates weakened 
magnetic correlations concomitant with "glassiness" and 
parasitic fM. In the following we focus on the type of 
samples showing the sharpest features in magnetization, 
where diffraction reveal sharp CO and magnetic reflec- 
tions. Above ~400 K the inverse susceptibility H/M (in- 
set Fig. [2k) follows a Curie- Weiss law with the effective 
moment /x e ff=5.51(9)/XB expected for Fe 2+ / 3+ and a neg- 
ative Weiss temperature of 6=— 307(9) K suggesting dom- 
inantly antiferromagnetic interactions, similar to YFe204 

Isothermal magnetization M(H) below Tjv (Fig.[3K/b) 
indicates a first-order metamagnetic transition, which be- 
comes strongly hysteretic for lower T. The hysteretic re- 
gion is indicated by the hatched area in Fig. QJi. The 
low-T saturation moment of the high- if phase (Fig. [2ji) , 
is similar to previous findings 16|, |21| , implying fM spin 
ordering. In contrast to this, the low-if phase near Tn 
seems to be AFM 25], with M oc if and no remanent 
moment. 

An AFM phase at 220 K and if = is inconsistent with 
the fM spin structure previously proposed [15j. The ab- 
sence of significant remanent M at 220 K could in princi- 
ple be explained by the formation of compensating fM do- 
mains. However, the drastic effect of if observed on sev- 
eral reflections with neutron diffraction, including a de- 
crease of (3, 3, integer) and increase of (|, |,halfmteger) 
reflections (Fig. [3^) and the emergence of additional in- 
tensity on structural reflections (Fig. |3t/d) show that 
the step in M(if) clearly corresponds to a coherent ef- 
fect, i.e., a genuine metamagnetic transition between two 
spin structures. The transition temperatures and fields 
from neutron scattering confirm the phase diagram from 
M(H,T), including the large hysteresis. 

The zero-field spin structure proposed in describes 
very well the (|, |, integer) reflections, but do not ac- 
count for newly observed magnetic reflections: With soft 
x-ray diffraction we observed a sharp reflection at (00 1) 
when the energy is tuned to the Fe L3 edge. According 
to previous work the polarization analysis (inset Fig. f2J:) 
suggests that this is purely magnetic, resulting, as ex- 
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FIG. 3. (Color online) (a,b) M(H) curves. The measurement 
direction is indicated by arrows. Virgin curves are measured 
after cooling in H — 0. (c,d) Integrated intensity of different 
reflections in neutron diffraction (D23) vs. H; compared with 
M{H), in (c) only the intensity change to 7q on structural 
reflections is illustrated. 



pected, from spins ||chex [27| . The similar T-dependence 
to the reflection (Fig. indicates that it origi- 

nates from the same spin structure, further supported by 
the suppression of the equivalent (102) + (00|) reflection 
at the metamagnetic transition (Fig. [3ji) . 

For the spin model in [l5| calculations show zero mag- 
netic intensity for these reflections. This model therefore 
has to be excluded. This spin structure 15j resulted from 
representation analysis based on the then only known 
i?3m crystallographic cell with no CO and a single Fe 
site, leading to a very small number of spin structures 
to be considered. To describe the new observed mag- 
netic reflections, we work within a 6x larger C2/m CO 
cell 14j, which corresponds to the magnetic cell for one 
domain according to all observed magnetic reflections. 

We take the most expansive approach by ignoring sym- 
metry and considering all 3 12 possible spin configurations 
of 12 Fe Ising spins (allowing for partial disorder) in the 
primitive cell |28[. Of these, ^15000 structures yield the 
same relative intensities at (|, |, integer) as the structure 
proposed in [l^|. To distinguish these structures, some 
broad-size restrictions for solutions can be made, based 
on the relative magnetic contribution of S + (00|) reflec- 

■| ,halfinteger) and struc- 
These restrictions show 



tions and an upper limit of 
tural S + (000) reflections [i_ 
that 7 symmetry-inequivalent spin structures can possi- 
bly be consistent with the observed magnetic diffraction 
in zero field. Refining these by fitting domain popula- 
tions and a Debye- Waller factor as in [ll|, but includ- 
ing (|, i,halfinteger), 6 models are rejected due to very 
large reduced \ 2 > the remaining solution is the AFM spin 
structure shown in Fig. [1}d. In contrast to the rejected 
structures, this solution is fully ordered and has a simple 
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FIG. 4. (color online), (a) Neutron diffraction (N5) pattern 
at 220 K along the (i|l) line in #=2.5T and in H=0. (b,c) 
Integrated intensity for scans along (||^) in both magnetic 
phases corrected at 220 K as described in the text. The solid 
red line represents the result from the model of [f5[ for the 
spin structures shown in Fig.f. The gray area indicate the 
magnet dark angle and the red arrow a reflection affected by 
a second grain. 



relationship to the high-TJ spin structure (see below) . In 
this structure the spins of each bilayer are fM aligned 
(tt4-)j hut the net moments of the bilayers are stacked 
antiferromagnetically, leading to the observed AFM be- 
havior. 

The moderate x 2 01 2.15 is due to systematically 
slightly higher intensities on (|, i,halfinteger) reflections 
(Fig. HJd). A cause for this could be magnetic contrast 
due to different Fe 2+ and Fe 3+ moments. For the CO 
proposed in 0, [l2| no significant improvement on refine- 
ment with Fe 2+ / 3+ magnetic contrast is observed for all 
possible spin structures, but a CO configuration previ- 
ously rejected due to charged bilayers [12| can further 
reduce x 2 to ~ 1 for the above optimal spin structure. 
Furthermore, a \ 2 ~1 can also be reached by considering 
a high-_ff phase contamination, albeit with a 15% phase 
fraction, which appears inconsistent with the remanent 
magnetization observed in Fig. |3]i. Given the similar 
effects of cross-contamination and CO, the Fe 2+ / 3+ dis- 
tribution cannot be conclusively [281 ] . 

The same approach was used for the high-i? phase [28[ . 
Comparing all with £ integer and halfintcger val- 

ues according to the diffraction pattern in 2.5 T (Fig.|3h) 
245 possibilities remain. The (T02) + (00§) reflection is 
strongly suppressed in the high-iJ phase (Fig.[3ji), reduc- 
ing the possibilities to 42. After comparing the magnetic 
contribution on different structural reflections (Fig. [3fc), 
only 18 solutions remain, corresponding to 3 symmetry 
inequivalent structures. Upon refinement, two are im- 
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FIG. 5. (Color online) Reciprocal space map of the scattered 
intensity (logarithmic scale) in the(hhl) plane by polarized 
neutrons in the spin flip channel (DNS) at different T. 



mediately rejected, the other (Fig. 0J;) is presented in 
Fig. [It. 

The fM solution is identical to the AFM solution, ex- 
cept that all Ising-spins in one Fe-0 bilayer flip their sign, 
leading to the overall 2 : 1 configuration of t and | spins 
consistent with the observed (Fig.^ji) net moment. This 
different stacking of bilayer net magnetization between 
the AFM (tit ■ • ■) and the fM (fit ■ • ■) phase shares 
similarity with the competing CO instabilities at higher 
T. Phase competition and metamagnetic transitions be- 
tween tit an d ttt-stacking of net moments are expected 
for layered magnets with very strong Ising-anisotropy [32I ] 
and has been observed in a few model systems at low T, 
e.g. FeCl 2 [IH. 

Intriguingly for LuFe2C>4, in contrast to expectations 
in simple model systems, the AFM-fM transition extrap- 
olates to H « for T — »■ Tjv as seen in Fig.[3}D; i.e. at Tjv 
and H = the two phases seem to be essentially degen- 
erate. The near-degeneracy of both charge and magnetic 
order is a hallmark of the importance of geometrical frus- 
tration in this system. The AFM/fM near-degeneracy in 
low H can lead to parts of the sample being trapped 
in fM after cooling through Tjf, particularly for samples 
with reduced T N (Fig.Qb). 

The particular differences between the two nearly- 
degenerate spin-structures suggest that the intra-bilayer 
correlations are more dominant than the inter-bilayer 
correlations. Just above TV we may therefore expect 
a random stacking of the net moment of still medium- 
range ordered bilayers, i.e. a 2D-order [29]. In contrast to 
the ferrimagnetic ordered bilayers of LuFe 2 04, for FeCl2 
the spins on triangular single layers are ferromagnetically 
coupled. For LuFe204, magnetic diffraction would result 
in strong diffuse scattering lines through {\\t) above Tn, 
still reasonably sharp in-plane, but featureless along £. 
This is indeed observed, visible at 280 K in Fig.[SJ Strong 
deviations from Curie- Weiss behavior up to ^400 K (in- 
set Fig. [5^) 0, , imply that these short-range corre- 
lations are significant in a wide T-range including Tco 
and may influence the establis hing of CO p}; provided 
there is a spin charge coupling [l7j . 

Although the AFM/fM metamagnetism presented here 
may not be resolvable in a majority of LuFe204 samples, 



the complex phase competition likely has ramifications 
for all specimens of this material. For example, if dis- 
order, e.g. due to oxygen off-stoichiometry, is added to 
the competing interactions, glassy freezing may be ex- 
pect to replace long-range spin order at TV, as observed 
in some samples [3, [23| ■ Disorder will disrupt most eas- 
ily the weak inter-bilayer correlations. It is thus natural 
to expect 3D spin order to be replaced by "spin-glass-like 
2D-ferrimagnetic order" , as reported from early neutron 
diffraction studies [2l|, and for the related YFe2 04-a; 
clearly linked to oxygen non-stoichiometry 26 1. 

In summary, we have elucidated the magnetic phase 
diagram of LuFe204 (Fig.[T]) close to T/v and determined 
a ferrimagnetic and an antiferromagnetic spin configura- 
tion, which are almost degenerate at TV and H = 0. This 
phase competition is not observed in classical metamag- 
netic materials, but is remarkably similar to competing 
CO instabilities at higher T. The near-degeneracy arises 
from geometrical frustration. Together with disorder it 
can lead for example to glassy freezing instead of long- 
range order, as observed in many samples, and competing 
magnetic fluctuations may influence the charge ordering 
and magnetoelectric coupling. 
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